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ABSTRACT 

A  post-coronal  mass  ejection  (CME)  current  sheet  (CS)  is  a  common  feature  developed  behind  an  erupting  flux  rope 
in  CME  models.  Observationally,  white  light  observations  have  recorded  many  occurrences  of  a  thin  ray  appearing 
behind  a  CME  eruption  that  closely  resembles  a  post-CME  CS  in  its  spatial  correspondence  and  morphology. 
UV  and  X-ray  observations  further  strengthen  this  interpretation  by  the  observations  of  high-temperature  emission 
at  locations  consistent  with  model  predictions.  The  next  question  then  becomes  whether  the  properties  inside  a 
post-CME  CS  predicted  by  a  model  agree  with  observed  properties.  In  this  work,  we  assume  that  the  post-CME  CS 
is  a  consequence  of  Petschek-like  reconnection  and  that  the  observed  ray-like  structure  is  bounded  by  a  pair  of  slow 
mode  shocks  developed  from  the  reconnection  site.  We  perform  time-dependent  ionization  calculations  and  model 
the  UV  line  emission.  We  find  that  such  a  model  is  consistent  with  SOHO/ UVCS  observations  of  the  post-CME 
CS.  The  change  of  Fexvm  emission  in  one  event  implies  an  inflow  speed  of  ~10  km  s-1  and  a  corresponding 
reconnection  rate  of  MA  ~  0.01 .  We  calculate  the  expected  X-ray  emission  for  comparison  with  X-ray  observations 
by  Hinode/X RT,  as  well  as  the  ionic  charge  states  as  would  be  measured  in  situ  at  1  AU.  We  find  that  the  predicted 
count  rate  for  Hinode /XRT  agrees  with  what  was  observed  in  a  post-CME  CS  on  2008  April  9,  and  the  predicted 
ionic  charge  states  are  consistent  with  high  ionization  states  commonly  measured  in  the  interplanetary  CMEs.  The 
model  results  depend  strongly  on  the  physical  parameters  in  the  ambient  corona,  namely  the  coronal  magnetic  field, 
the  electron  density,  and  temperature  during  the  CME  event.  It  is  crucial  to  obtain  these  ambient  coronal  parameters 
and  as  many  facets  of  the  CS  properties  as  possible  by  observational  means  so  that  the  post-CME  CS  models  can 
be  scrutinized  more  effectively. 
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1.  INTRODUCTION 

Most  models  for  coronal  mass  ejections  (CMEs),  regardless 
of  what  initiates  the  CME,  predict  a  current  sheet  (CS)  that 
develops  beneath  an  erupting  flux  rope  due  to  the  stretching 
of  the  overlying  coronal  field  (e.g.,  Lin  &  Forbes  2000;  Lynch 
et  al.  2004;  MacNeice  et  al.  2004;  Manchester  et  al.  2004). 
Magnetic  reconnection  in  the  CS  reduces  tension  restraining 
the  outgoing  flux  rope  and  at  the  same  time  produces  post- 
CME  loops  beneath  the  reconnection  point.  Outflows  along 
the  CS  and  a  temperature  higher  than  that  in  the  ambient 
corona  are  expected  inside  the  CS  as  magnetic  energy  is 
converted  to  kinetic  and  thermal  energy  due  to  reconnection. 
This  standard  flare-CME  picture  (Figure  1(a))  is  supported 
by  observations  such  as  loop-top  hard  X-ray  sources  (e.g., 
Masuda  et  al.  1994),  upward  growth  of  flare  loops  (e.g.,  Svestka 
1996),  separating  flare  ribbons,  and  hotter  post-flare  loops 
lying  higher  than  the  cooler  post-flare  loops  (e.g.,  van  Driel- 
Gesztelyi  et  al.  1997).  However,  few  observational  signatures 
of  the  CS  above  or  near  the  X-point  were  reported  until 
recently.  Sui  &  Holman  (2003)  and  Sui  et  al.  (2004)  reported 
signatures  of  a  CS  at  both  sides  of  the  reconnection  site  from 
X-ray  observations  by  the  Reuven  Ramaty  High  Energy  Solar 
Spectroscopy  Imager.  Savage  et  al.  (2010)  reported  a  post-CME 
CS  feature  observed  by  Hinode /XRT  with  downflows/upflows 
that  allow  the  X-point  to  be  located.  Innes  et  al.  (2003b)  and 
Wang  et  al.  (2007)  reported  high-speed  outflows  in  opposite 
directions  away  from  the  reconnection  site  based  on  Doppler 


shift  signatures  of  Fe  xix  and  Fe  xxi  lines  observed  by  Solar 
Ultraviolet  Measurements  of  Emitted  Radiation  (SUMER)  on 
the  Solar  and  Heliospheric  Observatory  (SOHO).  White  light 
(WL)  observations  of  CME  events  by  Solar  Maximum  Mission 
and  Large  Angle  Spectrometric  Coronagraph  (LAS CO)  on 
SOHO  have  recorded  many  occurrences  of  a  thin  ray  appearing 
behind  a  CME  eruption  (we  will  call  it  “WL  ray”  hereafter)  that 
closely  resembles  a  post-CME  CS  in  its  spatial  correspondence 
and  morphology  (e.g.,  Webb  et  al.  2003;  Ko  et  al.  2003;  Lin  et  al. 
2005).  Observations  by  Ultraviolet  Coronagraph  Spectrometer 
(UVCS)  on  SOHO  of  high- temperature  emission  (3-6  million 
degrees)  that  lies  along  such  WL  ray  in  between  the  erupting 
CME  and  post-CME  loops  provide  strong  support  for  the  CS 
interpretation,  as  opposed  to  being  just  a  usual  streamer  seen 
edge-on  (Ciaravella  et  al.  2002;  Ko  et  al.  2003;  Bemporad  et  al. 
2006;  Lee  et  al.  2006;  Ciaravella  &  Raymond  2008).  Figure  2 
shows  two  examples  of  such  events. 

Even  though  there  is  strong  observational  support  for  the 
observed  WL  ray  being  the  post-CME  CS,  that  interpretation 
remains  under  scrutiny.  One  major  cause  of  skepticism  arises 
from  the  observed  thickness  of  the  WL  ray,  as  well  as  the 
spatial  extent  of  the  high-temperature  emission  observed  by 
UVCS,  which  are  much  larger  than  what  most  reconnection 
models  predict.  The  observed  thickness  is  of  the  order  of 
104-105  km  (Lin  et  al.  2007,  2009;  Vrsnak  et  al.  2009)  which 
is  orders  of  magnitude  larger  than  theoretical  estimates  of  the 
thickness  for  the  diffusion  region  (DR)  and  Sweet-Parker  CS 
(see  Bemporad  2008;  Lin  et  al.  2009)  even  for  the  largest 
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(a)  (b)  (c) 

Figure  1.  (a)  Standard  flare-CME  model  depicting  a  CS  between  the  post-flare  loops  and  the  ejecting  flux  rope  (e.g.,  Lin  &  Forbes  2000).  (b)  Sketch  of  the  model  by 
Vrsnak  et  al.  (2009)  that  interprets  the  observed  WL  ray  (CS)  as  the  feature  bounded  by  a  pair  of  SMS  above  the  diffusion  region  (DR),  (c)  Sketch  of  the  “Fully-mixed 
model”  (adopted  from  Vrsnak  et  al.  2009). 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 


estimates  of  anomalous  resistivity  currently  available,  and 
taking  projection  effects  into  account.  Bemporad  (2008)  and 
Lin  et  al.  (2009)  discussed  several  reconnection  schemes  with 
turbulence  or  instability  (see  references  therein)  that  may 
produce  CS  thickness  comparable  to  what  was  observed,  such 
as  turbulence  and  stochastic/fractal/time-dependent  Petcheck- 
type  reconnection.  Model  calculations  have  shown  that  the 
effect  of  thermal  conduction  can  create  a  thermal  “halo”  around 
the  current  layer  (Yokoyama  &  Shibata  1997,  2001;  Seaton  & 
Forbes  2009).  This  would  contribute  to  the  extent  of  the  high- 
temperature  emission  observed  by  UVCS,  although  not  for  the 
WL  emission  (which  depends  mainly  on  the  electron  density). 
The  electric  potential  of  the  slow  mode  shocks  (SMSs)  might 
inhibit  transport  of  energy  into  the  upstream  flow,  reducing  this 
effect.  To  see  which  theories/models  are  more  consistent  with 
or  to  reject  the  CS  interpretation,  the  theories  or  models  for  the 
CS  must  predict  some  physical  quantities  that  can  be  tested  by 
these  observations. 

We  should  not  forget  that  there  is  another  piece  of  obser¬ 
vational  evidence  that  can  be  used  to  scrutinize  the  CS  inter¬ 
pretation  and  constrain  viable  theories /models.  As  mentioned 
above,  SOHO/ UVCS  has  observed  several  CME  events  ex¬ 
hibiting  high-temperature  emission  from  the  [Fexvui]  A974 
line  (formation  temperature  of  6  million  degrees)  at  heliocen¬ 
tric  heights  of  1.5-1. 7  Rq  that  lay  along  the  line  connecting 
the  erupting  CME  and  the  associated  post-CME  loops  (e.g.,  see 
examples  in  Figure  2).  Its  location  and  the  timing  of  its  bright¬ 
ening  following  the  eruption  strongly  favor  the  interpretation  of 
a  post-CME  CS.  It  is  then  important  to  go  one  step  beyond  mor¬ 
phology  and  find  out  if  a  given  CS  model  can  predict  physical 
quantities,  such  as  density  and  line  intensities,  that  agree  with 
the  observations.  And  if  so,  how  the  observed  quantities  would 
imply  the  physical  conditions  for  the  associated  reconnection 
within  the  context  of  a  given  model. 

It  is  important  to  note  that  the  analysis  of  emission  lines 
from  CSs  has  so  far  assumed  that  the  plasma  is  in  ioniza¬ 
tion  equilibrium,  but  this  can  be  a  poor  assumption.  For  in¬ 


stance,  the  ionization  time  scale  for  Fe  xvm  is  on  the  order  of 
1 010 / ne  s,  or  nearly  1000  s  for  the  densities  estimated  from 
UVCS  observations,  while  the  flow  time  in  the  CS  should  be  on 
the  order  of  a  few  tenths  of  a  solar  radius  divided  by  an  outflow 
speed  on  the  order  of  1000  km  s-1,  which  can  be  several  times 
smaller.  It  is  therefore  important  to  consider  time-dependent 
ionization  when  interpreting  emission  line  intensities. 

Vrsnak  et  al.  (2009)  proposed  a  working  hypothesis  for  the 
observed  post-CME  CS  (i.e.,  WL  ray),  and  provided  the  first 
attempt  to  calculate  the  electron  density  within  the  CS,  and 
compare  with  the  observations.  This  model  assumes  a  steady- 
state  Petschek-like  reconnection  scheme  (Petschek  1964).  The 
reconnection  takes  place  in  the  DR  and  a  pair  of  standing  SMSs 
develop  and  extend  out  from  the  DR.  At  the  SMS  crossing, 
the  inflowing  coronal  material  is  compressed  and  heated  in  a 
way  that  is  mainly  determined  by  the  external  plasma  /3,  and 
the  velocity  of  the  outflowing  jet  within  the  region  bounded 
by  the  SMSs  is  approximately  equal  to  the  external  Alfven 
speed  (Aurass  et  al.  2002).  They  found  that  the  observed 
electron  density  and  morphology  of  the  post-CME  CS  (including 
thickness)  can  be  successfully  explained  by  their  model.  In  this 
paper,  we  adopt  their  model  and  calculate,  for  the  first  time, 
the  expected  UV  and  X-ray  emission  inside  the  CS.  Strictly 
speaking,  within  the  context  of  this  model,  the  so-called  post- 
CME  CS  seen  as  a  ray-like  structure  in  WL  observations  is 
actually  the  region  bounded  by  the  SMSs.  This  region  is  not  a 
“current  sheet”  per  se,  and  the  electric  currents  are  concentrated 
only  in  the  DR  and  SMSs  (Vrsnak  et  al.  2009).  Note  that, 
however,  we  will  use  “post-CME  CS”,  “CS,”  or  “WL  ray” 
interchangeably  throughout  this  paper  to  stand  for  this  region. 

Section  2  describes  the  models  in  detail.  We  adopt  two 
approaches.  One  is  the  “Fully-mixed”  model  which  is  the  direct 
adaptation  of  the  model  laid  out  in  Vrsnak  et  al.  (2009).  The 
other  is  the  “Streamline  model”  which  assumes  that  the  outflow 
jets  out  of  the  SMS  crossings  at  different  locations  do  not  mix 
with  each  other.  We  describe  the  methods  for  calculating  the 
electron  density,  electron  temperature,  and  time-dependent  ionic 
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Figure  2.  Upper  panels:  the  CME/CS  event  on  2002  January  8  (see  Ko  et  al.  2003).  Lower  panels:  the  CME/CS  event  on  2003  November  4  (see  Ciaravella  & 
Raymond  2008).  Images  are  from  SOHO/ LASCO  (coronal  images  in  red-orange),  SOHO/ EIT  A.  1 95  (solar  images  in  green),  SOHO/ UVCS  (slit  images  marked  in 
the  two  right  panels),  SOHO /CDS  (coronal  loop  image  in  gray  in  the  upper  right  panel),  and  MLSO/MK4  (coronal  image  in  blue  in  the  upper  right  panel). 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 


fractions  within  the  CS  for  both  models.  Section  3  presents  the 
predicted  UV  and  X-ray  emission,  and  the  expected  frozen-in 
charge  states  as  would  be  measured  in  situ  at  1  AU.  We  discuss 
the  results  in  Section  4. 


2.  MODELING  THE  POST-CME  CURRENT  SHEET 

Figure  1(b)  illustrates  the  idea  of  the  model  of  Vrsnak  et  al. 
(2009)  that  is  based  on  the  steady-state  Petschek  reconnection 
scheme.  As  the  eruption  stretches  the  oppositely  directed  coro¬ 
nal  field  beneath  the  flux  rope,  reconnection  takes  place  at  the 
DR  under  suitable  conditions.  A  pair  of  standing  SMSs  then  de¬ 
velops  out  of  the  DR  both  above  and  below.  In  this  paper,  we  will 
concentrate  on  the  region  bounded  by  the  pair  of  SMSs  above 
the  DR  in  this  vertical  CS  configuration.  Coronal  material  that 
crosses  the  SMS  is  compressed,  heated,  accelerated,  and  forms 
an  outflow  jet  above  the  DR.  The  extent  of  compression  and 
heating  depends  on  the  ambient  (i.e.,  coronal)  plasma  /3  at  the 
location  of  the  SMS  crossing  as  in  Aurass  et  al.  (2002): 

ni  =  5(l  +  fl 

m  2  +  5/3  { 


^  =  i  +  A, 
Ti  5/3’ 


(2) 


where  n\,T\  and  are  the  electron  density  and  electron 

temperature  in  the  ambient  corona  and  after  the  SMS  crossing, 
respectively.  /3  is  the  ratio  of  gas  to  magnetic  pressure  in  the 
ambient  corona  for  fully  ionized  plasma  (Mann  et  al.  1999): 


l.92n\kT\ 
P  (2&/&T)  ’ 


(3) 


where  BCOT  is  the  ambient  coronal  magnetic  field  at  the  SMS 
crossing.  The  outflow  speed  out  of  the  SMS  crossing  is  ap¬ 
proximated  to  be  the  Alfven  speed  for  the  inflowing  material: 

VA  =  /y|  =  .  (4) 

^/47r  x  \.Yln\mp 

To  calculate  the  physical  properties  inside  this  post-CME  CS, 
one  needs  to  take  into  account  all  outflows  from  the  SMS 
crossings  along  the  SMS  (Figure  1(b)).  At  this  point,  we 
will  investigate  two  models.  One  is  the  “Fully-mixed”  model 
(“Model  1”,  Figure  1(c))  in  which  plasma  that  just  comes  inside 


628 


KO  ET  AL. 


Vol.  722 


the  CS  through  SMS  mixes  thoroughly  with  material  that  comes 
from  below,  as  in  Vrsnak  et  al.  (2009).  This  is  the  “fluid” 
approach  and  the  electron  conductivity  is  large  across  the  CS 
so  that  thermal  equilibrium  is  achieved  between  the  incoming 
and  outflow  material  at  a  given  height.  The  other  is  an  extreme 
opposite  of  the  first  model,  the  “Streamline”  model  (“Model  2”, 
cf.  Figure  1(b)),  in  which  plasma  that  just  comes  inside  the  CS 
through  SMS  forms  its  own  “streamline”  and  does  not  mix  with 
the  material  that  comes  in  at  other  SMS  crossings.  This  would 
be  the  case  if  certain  conditions  inside  the  CS,  e.g.,  turbulence, 
prevent  these  flows  (i.e.,  streamlines)  from  efficiently  changing 
the  thermodynamical  properties  of  each  other.  We  discuss  the 
two  models  in  more  detail  below. 

Both  models  adopt  the  same  ambient  coronal  conditions.  As 
shown  above,  the  properties  inside  the  CS  depend  on  the  ambient 
plasma  ft,  which  in  turn  depends  on  the  electron  density,  electron 
temperature,  and  magnetic  field  profiles  in  the  ambient  corona. 
For  the  coronal  magnetic  field,  we  adopt  the  empirical  model 
by  Dulk  &  McFean  (1978): 

Bcor(r)  =  0.5(r  -  l)"1'5  G,  (5) 

where  r  is  the  heliocentric  distance  from  the  Sun  in  Rq. 
Strictly  speaking,  this  formula  is  applicable  for  approximately 
1.02  ^  r  ^  10  Rq  but  we  will  use  this  formula  up  to  20  Rq 
(see  also  Vrsnak  et  al.  2002),  the  upper  boundary  of  our 
calculations.  The  accuracy  of  this  formula  between  10  and 
20  R0  is  n°t  a  concern  for  our  study  since,  as  we  will  see 
below,  (1)  the  emission  calculated  here  only  has  observational 
data  available  at  locations  much  lower  than  10  Rq  (because  the 
emission  drops  with  n2e)  and  (2)  the  change  of  time-dependent 
ionic  fractions  occurs  mostly  below  10  Rq.  This  formula  gives 
a  value  of  1.6  x  10-4  G  (16  nT)  at  1  AU  which  is  also 
reasonable  (e.g.,  see  ACE  E3  summary  plots  at  the  ACE  Science 
Center,  http://www.srl.caltech.edu/ACE/ASC/index.html).  To 
investigate  the  dependence  on  different  coronal  electron  density 
and  temperature  profiles,  we  use  two  nCOY  and  two  Tcor  profiles. 
For  the  coronal  electron  density  profiles,  one  (denoted  as  “Nl”) 
is  a  “hybrid”  profile  which,  for  r  ^  5.66  Rq,  is  obtained  for  a 
streamer  during  the  SPARTAN  201-1  mission  (Guhathakurta  & 
Fisher  1995).  Beyond  5.66  RQ ,  we  linearly  interpolate  between 
the  measured  values  at  5.66  Rq  and  a  value  of  211.2  cm-3  at 
43  Rq  (0.2  AU)  which  is  based  on  Mann  et  al.  (1999): 

ncor(r)  =  6.53  (  — )  cm~3  (6) 

for  0.2  AU  ^  r  ^  5  AU.  The  second  ne(r)  (denoted  as  “N2”)  is 
adopted  from  LeBlanc  et  al.  (1998)  which  is  of  an  analytic  form 
of 

nmr(r)  =  3.3  x  105r-2  +4.1  x  10 6r~4  +  8.0  x  lOV^cm-3. 

(7) 


Note  that  both  profiles  are  based  on  observations.  For  the  coronal 
electron  temperature  profiles,  we  assume  a  general  form  (in  K): 


|  To,  for  1  Rq  ^  r  ^  2  Rq 

{:To/(0.5r)0-92  for  r  >  2  Re. 


(8) 


The  index  0.92  is  based  on  Totten  et  al.  (1995)  for  a  spherically 
expanding  wind.  For  Tq,  we  will  compare  results  for  Tq  = 
106  K  (71)  and  2  x  106  K  ( T2 ).  The  solar  wind  speed  profile  is 
then  taken  as 

o  1  /  215\2 

vsw(r)  =  2  x  108 - —  -  cm  s  1  (9) 

nCov(r)  \  r  ) 


based  on  the  continuity  equation  nvr 2  =constant  with  nv  — 
2  x  108  cm-2  s-1  at  1  AU. 

Figure  3  plots  the  adopted  Bc or,  nc or,  and  Tcor  profiles,  and 
the  corresponding  plasma  Also  plotted  (“lMK_Parker”)  are 
those  adopted  by  Vrsnak  et  al.  (2009)  using  a  Parker  wind 
model  (Mann  et  al.  1999)  with  1  MK  isothermal  corona,  and 
the  same  Bcor(r)  (Dulk  &  McFean  1978).  Note  that  none  of 
the  actual  coronal  parameters  in  the  immediate  vicinity  of  the 
SMSs  are  known  for  certain,  and  the  coronal  conditions  are 
likely  to  vary  from  one  event  to  another.  Furthermore,  Vrsnak 
et  al.  (2009)  found  that  the  coronal  density  is  depleted  from  the 
pre-CME  corona  in  the  vicinity  of  the  WF  ray.  In  any  case,  we 
believe  that  these  coronal  profiles  adopted  here  are  reasonable 
and  can  be  taken  as  generic  choices.  Nonetheless,  it  is  necessary 
to  investigate  the  effect  of  different  coronal  parameters  on  the 
calculated  CS  properties.  Therefore,  we  calculate  the  model 
results  for  three  cases:  (1)  Nl+Tl  (denoted  as  “N1T1”),  (2) 
N1+T2  (denoted  as  “N1T2”),  and  (3)  N2+T2  (denoted  as 
“N2T2”)  while  keeping  the  same  BCOY(r).  Cases  “N1T1”  and 
“N1T2”  are  used  to  investigate  effects  from  different  coronal 
temperature  profiles.  Cases  “N1T2”  and  “N2T2”  are  used  to 
investigate  effects  from  different  coronal  density  profiles.  We 
take  N1T2  as  the  reference  case  for  comparison  among  different 
cases  and  models.  The  lower  right  panel  of  Figure  3  plots  the 
plasma  which  determines  the  jump  condition  across  the  SMS 
(Equations  (1)  and  (2)),  for  these  three  cases.  The  in  Vrsnak 
et  al.  (2009,  “lMK_Parker”)  is  also  plotted  for  comparison. 

2.7.  Model  1:  The  Fully -mixed  Model 

Model  1  (Figure  1(c))  follows  the  same  thermodynamic 
calculations  as  Vrsnak  et  al.  (2009),  except  that  we  calculate  the 
CS  properties  with  three  cases  of  coronal  parameters  described 
above.  Readers  are  referred  to  the  appendix  in  Vrsnak  et  al. 
(2009)  for  a  complete  description  of  the  model.  We  only  briefly 
describe  the  concept  here.  For  a  given  cell  (A R  in  Figure  1(c)) 
at  a  given  height  Rt  above  the  DR,  plasma  comes  into  the 
CS  through  the  SMS  with  density  n 2  (Equation  (1))  into  a 
given  volume  (Vsms)  that  is  determined  by  the  ambient  coronal 
parameters.  There  is  also  material  coming  from  the  cell  below 
into  a  volume  (Vthm)  determined  by  spherical  expansion.  In 
this  model,  the  outflow  (i.e.,  exhaust  flow)  speed  is  taken  as 
the  local  value  of  va  +  iw  (for  details  see  Vrsnak  et  al.  2009) 
but  usw  is  almost  always  negligible  in  the  considered  height 
range.  Under  the  assumption  that  the  “incoming”  and  “through” 
material  thoroughly  mix  with  each  other,  the  electron  density 
flowing  out  of  this  cell  into  the  cell  above  it  can  be  calculated 
from  mass  conservation,  along  with  the  CS  geometry.  The 
electron  temperature  in  the  cell  is  calculated  to  be  the  average  of 
the  temperature  after  the  SMS  crossing  (i.e.,  T2,  Equation  (2)) 
and  that  coming  from  the  cell  below  under  adiabatic  cooling, 
weighted  by  the  mass  inside  Vsms  and  Vthm,  respectively.  The 
density  and  temperature  profiles  along  the  CS  (above  the  DR) 
can  then  be  calculated  iteratively,  and  they  depend  on  the 
ambient  coronal  parameters  and  the  location  of  the  DR.  Since 
the  magnetic  field  in  this  configuration  is  perpendicular  to  the 
exhaust  flow,  the  thermal  conduction  in  the  flow  direction  will 
be  inhibited.  Note  that  the  inflow  speed  V{n,  which  is  a  factor  that 
governs  the  CS  geometry,  is  calculated  from  r  vmBcov  =  constant 
under  the  steady-state  assumption  and  with  V{n  =  0.01  va  at 
the  DR. 

Figure  4  plots  the  resulting  ne(r ),  Te(r),  and  vollt(r)  for  the 
three  cases  of  the  coronal  ne/ Te  profiles,  and  for  DR  at  1.1  and 
1.5  Rq  for  each  case.  The  profiles  from  Vrsnak  et  al.  (2009, 
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Figure  3.  Profiles  for  coronal  magnetic  field  strength  (upper  left),  electron  density  (upper  right),  electron  temperature  (lower  left),  and  plasma  beta  (lower  right)  adopted 
for  the  study.  Also  plotted  is  the  electron  density  and  plasma  beta  profiles  for  the  “lMK_Parker”  model  adopted  in  Vrsnak  et  al.  (2009).  The  electron  temperature 
profile  for  the  lMK_Parker  model  is  isothermal  at  Te(r)  =  1  x  106  K. 


“lMK_Parker”)  and  the  electron  densities  derived  from  UVCS 
and  LASCO  data  are  also  plotted  for  comparison.  The  profiles 
with  DR  heights  in  between  1.1  and  1.5  RQ  lie  between  the 
two  curves  for  each  case,  respectively.  Note  that  models  N1T1 
and  N1T2  have  almost  the  same  profiles  in  the  CS  because 
when  beta  is  small  as  in  these  cases,  the  jump  condition  for  ne 
is  about  the  same  (Equation  (1)).  The  factor  of  two  difference 
in  is  compensated  by  the  factor  of  two  in  T\  (Equation  (2) 
implies  T2/T\  ~  for  small  /3),  leading  to  almost  the  same 
temperature  T2  in  the  CS.  The  electron  density  in  the  CS  is 
higher  than  that  in  the  ambient  corona  at  all  heights  (comparing 
Figures  3  and  4;  see  also  Equation  (1));  therefore,  the  CS 
structure  can  stand  out  against  the  ambient  corona  as  observed. 

In  order  to  calculate  the  emission  from  the  CS,  we  need  to  first 
calculate  the  ion  charge  states  (i.e.,  ionic  fractions)  of  the  ions 
that  contribute  to  the  emission.  With  large  outflows  in  the  CS,  the 
ion  charge  states  do  not  always  maintain  ionization  equilibrium 
when  the  electron  density  is  low  enough.  The  evolution  of  the 
ion  charge  states  with  the  flow  depends  on  the  electron  density 
and  temperature  profiles  as  well  as  the  ion  outflow  profiles.  For 
a  given  element,  the  evolution  of  ion  charge  states  with  the  flow 
is  expressed  as  (under  steady- state  assumption): 


^ e  —  q  —  \  C i,qyq )  +  ij^rr,q+\  +  Rdr,q+ 1)  3^+1 

fie  (j^rr,q  +  Rdr,q )  37#>  (10) 


where  r  is  the  heliocentric  height  of  the  fluid  element  along  the 
CS,  yq  is  the  ionic  fraction  of  the  charge  state  q.  C^q ,  Rrr,q ,  Rdr,q 
(which  mainly  depend  on  Te )  are  the  rates  for  electron  impact 
ionization  (including  auto-ionization),  radiative  recombination, 


and  dielectronic  recombination,  respectively,  out  of  the  charge 
state  q.  The  ionization  and  recombination  rates  are  calculated 
using  the  most  recent  compilation  by  Bryans  et  al.  (2006,  2009, 
and  references  within).  Here  we  assume  that  all  ions  have  the 
same  outflow  speed.  For  each  element,  Equation  (10)  for  all 
charge  states  (i.e.,  Z+l  equations  for  an  element  with  atomic 
number  Z)  is  solved  simultaneously.  We  calculate  the  ionic 
fractions  for  13  most  abundant  elements:  H,  He,  C,  N,  O,  Ne, 
Mg,  Si,  S,  Ar,  Ca,  Fe,  and  Ni.  Note  that  the  ions  that  just 
cross  the  SMS  should  still  carry  the  charge  state  distribution 
at  the  ambient  coronal  temperature  (even  though  the  electron 
temperature  jumps  to  T2  after  the  crossing),  and  those  calculated 
by  Equation  (10)  are  for  ions  coming  from  the  cell  below  (cf. 
Figure  1(c))  following  the  flow.  Therefore,  the  charge  state 
distribution  in  a  given  cell  would  be  the  average  of  the  two 
regions  weighted  by  the  mass  inside  Vsms  and  Vthm»  respectively, 
for  that  given  cell. 

Figure  5  plots  the  resulting  yq(r)  for  case  N1T2  with  DR  at 
1 .2  Rq  for  a  few  ions,  along  with  the  values  in  the  case  of  ioniza¬ 
tion  equilibrium.  We  can  see  that  these  ion  charge  states  are  far 
from  ionization  equilibrium  above  a  certain  height  (different  for 
different  ion  species),  and  the  ionic  fractions  “freeze-in”  as  they 
do  in  the  solar  wind.  Therefore,  to  predict  emissions  inside  the 
CS,  it  is  important  to  take  this  non-equilibrium  condition  into 
account,  and  it  is  not  necessarily  valid  to  calculate  the  emission 
based  on  ionization  equilibrium  (i.e.,  only  based  on  the  elec¬ 
tron  temperature).  Similarly,  the  electron  temperature  derived 
from  line  ratios  assuming  ionization  equilibrium  may  not  be  the 
actual  electron  temperature. 

Figure  6  compares  yq(r)  for  the  three  cases  at  DR  height 
of  1.2  Rq  for  a  few  ions.  We  can  see  that  there  are  significant 
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(Figure  4),  the  different  ambient  coronal  temperatures  result  in 
different  initial  charge  states  at  the  location  of  the  SMS  crossing, 
thus  different  evolution  of  the  ionic  fractions  as  the  ions  flow 
along  the  CS.  Therefore,  the  observations  can,  in  principle,  be 
used  to  constrain  these  input  model  parameters.  Note  that  many 
ions  seem  to  start  freezing-in  at  higher  heights  than  what  are 
usually  expected  for  the  normal  solar  wind  (e.g.,  Biirgi  &  Geiss 
1986;  Ko  et  al.  1997),  probably  due  to  much  higher  density 
inside  the  CS  than  in  the  solar  wind. 

2.2.  Model  2:  The  Streamline  Model 

Another  model  we  explored  assumes  that  the  plasma  that  just 
comes  inside  the  CS  through  SMS  forms  its  own  “streamlines” 
and  does  not  mix  with  the  material  that  comes  in  at  other  SMS 
crossings.  The  properties  inside  the  SMS  at  a  given  location 
would  then  be  an  average  quantity  from  these  streamlines 
weighted  by  the  volume  each  of  these  streamlines  occupies  at 
that  location  (see  below).  While  the  long  collisional  mean  free 
path  suggests  the  idea  of  mixing  as  in  Model  1 ,  it  is  not  clear  if  the 
actual  magnetic  field  configuration  inside  the  SMS  is  the  same  as 
that  in  the  Petchek’s  model.  Also,  MHD  turbulence  may  inhibit 
the  mixing  of  the  plasma  and  thermal  equilibration  (turbulent 
speeds  of  ^60  km  s-1  have  been  found  in  the  post-CME  CSs  by 
Bemporad  2008).  Therefore,  we  also  want  to  examine  how  such 
a  “Streamline  model”,  which  can  be  regarded  as  an  extreme 
opposite  case  of  Model  1,  compares  with  the  observations. 

We  use  the  same  three  coronal  models  (i.e.,  N1T1,  N1T2, 
N2T2),  and  the  same  compression/heating/acceleration  at  the 
SMS  crossing  (Equations  (1),  (2),  and  (4)).  As  the  plasma  flows 
in  each  streamline  after  the  SMS  crossing  (cf.  Figure  1(b)),  the 
electron  density  decreases  due  to  spherical  expansion: 


Figure  4.  ne(r),  Te(r),  and  vout(r)  for  the  three  cases  of  the  coronal  ne/Te 
profiles,  and  for  DR  at  1.1  (in  black)  and  1.5  Rq  (in  red)  for  each  case. 
The  profiles  from  Vrsnak  et  al.  (2009,  “lMK_Parker”;  dash-triple-dot)  and 
the  electron  densities  derived  from  UVCS  and  LASCO  data  (blue  bars  and 
pluses)  are  also  plotted  for  comparison.  UVCS  data  are  from  the  1998  March  23 
event  (Ciaravella  et  al.  2002)  at  1.50  Rq,  2002  November  26  event  (Bemporad 
et  al.  2006)  at  1.61  Rq,  and  the  2003  November  4  event  (Ciaravella  &  Raymond 
2008)  at  1.66  R0-  LASCO  data  are  from  the  2002  January  8  event  (Ko  et  al. 
2003)  at  3.0  and  4.4  Rq.  Note  that  since,  with  the  same  coronal  field,  uout  only 
depends  on  the  ambient  density,  it  is  the  same  for  cases  N1T1  and  N1T2. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 

differences  between  models  of  different  coronal  profiles  that 
would  result  in  different  emission  intensities  at  various  heights, 
as  well  as  different  frozen-in  charge  states.  Note  that  even  though 
models  N1T1  and  N1T2  have  almost  the  same  ne ,  Te ,  and  vout 


r 2 

ne(r)  =  «e,sMS-^,  (11) 

where  ne, sms  (same  as  n2  in  Equation  (1))  is  the  electron 
density  at  the  SMS  crossing  at  height  tsms-  The  evolution  of 
the  electron  temperature  is  governed  by  adiabatic  cooling  and 
radiative  cooling: 

dTJf)  _  4  7; 

dr  3  r  3(1.92 neku{r)Y 

where  A  is  the  radiative  cooling  rate.  Because  the  ions  are 
most  likely  not  in  ionization  equilibrium  along  the  flow  (cf. 
Equation  (10)  and  Figure  5),  we  use  the  non-equilibrium  ionic 
fractions  (Equation  (10))  to  calculate  the  radiative  cooling 
rate.  The  radiative  cooling  includes  bound-bound,  bound-free, 
and  free-free  processes  using  routines  provided  by  CHIANTI 
atomic  database  version  6.0  (Dere  et  al.  1997,  2009)  but  with 
ionic  fractions  calculated  here  replacing  those  in  ionization 
equilibrium.  We  find  that  the  radiative  cooling  is  negligible 
compared  to  adiabatic  cooling  in  almost  all  cases  except  at  the 
few  lowest  heights  where  the  density  is  high.  Therefore  it  can  be 
neglected  in  general.  The  flow  speed  in  each  streamline  is  equal 
to  the  coronal  Alfven  speed  at  the  SMS  crossing  (Equation  (4)). 
The  evolution  of  the  ionic  fractions  in  each  streamline  thus  can 
be  calculated  according  to  Equation  (10).  Note  that  we  neglect 
heat  conduction  in  Equation  (12)  assuming  that  MHD  turbulence 
inside  the  CS  suppresses  the  thermal  conduction.  Long  CSs  are 
expected  to  undergo  tearing  instability,  which  can  create  a  large 
number  of  plasmoids  at  different  scales  (e.g.,  Shibata  &  Tanuma 
2001;  Barta  et  al.  2010,  and  references  therein).  Due  to  the 
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Figure  5.  Ionic  fraction  profiles  for  DR  height  at  1.2  Rq  for  the  Fully-mixed  model  with  coronal  model  N1T2,  compared  with  those  in  ionization  equilibrium  at  the 
local  electron  temperature  in  the  CS  (dotted,  dashed,  dash-dotted  in  the  order  of  ascending  charge  state). 


poloidal  field  of  the  plasmoids,  they  are  thermally  “insulated” 
in  such  medium  that  is  filled  with  large  number  of  stochastically 
moving  plasmoids,  thus  reducing  the  thermal  conduction.  There 
will  be  thermal  conduction  at  some  level  in  spite  of  the  effects  of 
turbulence,  but  we  have  chosen  the  extreme  case  to  investigate 
the  range  of  possible  results. 

Once  the  electron  density,  temperature,  and  ion  outflow 
profiles  in  each  streamline  are  known,  the  ionic  fraction  profiles 
in  the  streamline  can  be  calculated  according  to  Equation  (10). 
The  properties  inside  the  CS  at  a  certain  location,  in  the  context 
of  this  model,  are  therefore  average  quantities  taking  into 
account  all  “streamlines”  that  form  between  the  DR  and  that 


given  location,  assuming  that  the  streamline  structure  inside  the 
CS  cannot  be  resolved.  Under  spherical  geometry,  the  volume 
at  height  r  occupied  by  a  streamline  that  entered  the  CS  (i.e., 
crossing  the  SMS)  at  height  rSMS  is  proportional  to  r/rSMS-  The 
average  quantity  of  a  parameter  X  observed  at  height  r  is  thus 


X(r)  = 


rSMS) 

KZl=rmF(r,rSMS) 


(13) 


where  tSms  and  rDR  are  the  height  of  the  SMS  crossing  and 
DR,  respectively.  X$ms  is  the  physical  quantity  at  height  r  in 
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Figure  6.  Ionic  fraction  profiles  for  DR  height  at  1.2  Rq  for  the  Fully-mixed  model.  Plotted  are  for  three  coronal  models,  N1T1  (green  dash),  N1T2  (black  solid),  and 
N2T2  (red  dash  dot). 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 


the  streamline  that  enters  the  SMS  at  tsms>  and  F(r,  tsms)  is 
the  “contribution  factor”  from  the  streamline  forming  at  rSMS- 
F(r,  rs ms)  is  different  for  different  physical  parameters.  For  the 
electron  density,  F(r,r sms)  =  D^sms-  For  the  ionic  fraction 
yq,  F(r,  rSMs)  =  rne(r)/rsMS  (ion  density  is  proportional  to 
yq  *  ne). 

Figure  7  shows  an  example  of  how  different  streamlines 
contribute  to  an  ionic  fraction  observed  at  height  r  =  1.7  Rq 
for  the  case  of  rDR  =  1.2  RQ.  All  three  coronal  models  are 
plotted  for  comparison.  Also  plotted  are  those  for  the  electron 
density.  We  choose  those  ions  that  emit  high-temperature 
lines  observable  by  UVCS  in  the  CS  at  this  height  (i.e., 
Fe+17([Fe  xviii]  A974),  Si+11(Sixn  A499),  Ca+13([Caxiv]  A943) 
(e.g.,  Ciaravella  et  al.  2002;  Ko  et  al.  2003;  Bemporad  et  al. 
2006).  One  can  see  that  the  streamlines  (i.e.,  rSMs’s)  that  make 
a  major  contribution  to  the  average  ionic  fraction  are  different 
for  different  ions  and  different  coronal  models,  but  most  of 
the  contribution  is  from  SMS  crossings  at  low  heights.  Similar 
examination  for  an  observed  height  of  1  AU  indicates  that  the 
major  contribution  for  most  abundant  ions  at  1  AU  is  from 
streamlines  that  enter  the  SMS  below  2  RQ.  For  the  electron 
density,  it  is  the  same  for  all  three  models  because  ne(r)  evolves 
in  the  same  way  with  r  (Equation  (11)). 


Figure  8  plots  the  average  ionic  fractions  (Equation  (13)) 
for  selected  ions  as  a  function  of  height  for  the  three  coronal 
models.  Comparing  with  the  Fully-mixed  model,  we  find  that 
the  results  of  the  two  models  are  similar  within  an  order 
of  magnitude.  Figure  9  plots  the  average  electron  density 
compared  with  the  data  (cf.  Figure  4).  Figure  9  shows  that, 
similar  to  the  Fully-mixed  model,  models  N1T1  and  N1T2  with 
reasonable  DR  heights  predict  electron  densities  that  agree  with 
the  observations,  while  model  N2T2  seems  to  underestimate  the 
electron  density. 

3.  PREDICTED  OBSERVABFES  FROM  MODEF 
C  AFCUF  ATION  S 

The  main  purpose  of  this  work  is  to  test  post-CME  CS  models 
and  see  if  model  predictions  agree  with  observations.  Figures  4 
and  9  show  that  the  predicted  electron  density  profiles  within 
the  CS  overall  agree  with  the  data  for  both  model  approaches 
(see  also  Figure  7  of  Vrsnak  et  al.  2009).  The  two  models 
predict  similar  electron  densities  at  low  heights,  but  the  electron 
density  profile  falls  faster  with  height  for  the  Fully-mixed  model 
at  r  >  4  Rq.  Therefore,  CS  observations  over  an  extended 
range  of  heights  may  be  able  to  distinguish  such  different  model 
approaches  and  indicate  which  one  represents  better  the  actual 
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Figure  7.  Contribution  to  the  total  ionic  fraction  observed  at  height  1.7  Rq  from  streamlines  entering  the  SMS  at  heights  from  1.2  Rq  (DR  height)  to  1.7  Rq.  Plotted 
are  for  three  ions  and  the  electron  density  for  three  coronal  models,  N1T1  (green  dash),  N1T2  (black  solid),  N2T2  (red  dash  dot).  The  SMS  height  step  is  0.01  Rq. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 


conditions  inside  the  CS.  It  is  also  possible,  in  the  context  of 
either  model,  to  constrain  the  coronal  parameters  outside  of 
the  CS. 

3.1.  UV  Line  Emission 

Another  predicted  quantity  is  the  line  emission  (see 
Section  1).  In  coronal  conditions,  the  electron  density  is  low 
enough  that  most  emission  lines  are  produced  by  electron  col- 
lisional  excitation  followed  by  spontaneous  emission.  For  a 
plasma  with  electron  temperature  Te,  the  line  emission  is  thus 

/line=-L^  f  G(Te)dEM{Te)  photon  cm-2  s-1  sr-1, 

4 n  nH  J 

(14) 

where  ne\/nH  is  the  elemental  abundance  relative  to  hydrogen 
(i.e.,  absolute  abundance;  Grevesse  et  al.  2007).  G(Te)  is  the 
contribution  function  which  is  defined  as 

G(Te)  =  —Blmcqrmc(Te).  (15) 

nQ  i 

where  n[on/nQ\  is  the  ionic  fraction  calculated  from 
Equation  (10).  B\\ne  is  the  branching  ratio  for  the  line  transi¬ 
tion,  and  <7iine  ( )  is  the  electron  excitation  rate  which  is  a  func¬ 
tion  of  only  the  electron  temperature  in  the  low-density  limit. 
dEM(Te )  =  d(nennL)  is  the  emission  measure  (in  cm-5)  at  a 
given  electron  temperature  with  line-of-sight  (LOS)  depth  L,  and 
nH  =  ne/ 1.2  for  fully  ionized  plasma.  In  our  models,  the  ionic 


fractions  are  not  in  ionization  equilibrium  (e.g.,  see  Figure  5). 
Therefore,  we  calculate  G  from  the  CHIANTI  atomic  database 
version  6.0  (Dere  et  al.  1997,  2009)  but  with  the  ionic  fractions 
calculated  here  replacing  those  in  ionization  equilibrium  (e.g., 
Figures  6  and  8). 

For  the  Fully-mixed  model,  the  line  emission  can  be  cal¬ 
culated  directly  from  Equation  (14)  with  the  electron  den¬ 
sity,  temperature,  and  ionic  fractions  calculated  in  the  model. 
For  the  Streamline  model,  one  needs  to  calculate  the  average 
quantity  as  in  Equation  (13)  with  2fSMs(r)  =  7iine,SMs(r)  and 
F(r,  rs ms)  =  i* / ^sms-  There  are  two  assumptions  needed.  One 
is  the  elemental  abundance,  and  the  other  is  the  LOS  depth.  The 
first  ionization  potential  (FIP)  effect  is  commonly  observed  in 
coronal  plasmas.  The  FIP  bias,  defined  as  the  abundance  ratio  of 
a  low-FIP  element  (FIP  <10  eV)  to  a  high-FIP  element  (FIP  > 
10  eV)  relative  to  its  photospheric  ratio,  is  generally  in  the 
range  of  3-4  in  active  regions  and  streamers  (e.g.,  see  review  by 
Raymond  et  al.  2001,  and  references  within).  Since  fast  CMEs 
with  such  WL  ray/CS  structure  almost  always  occur  from  active 
regions,  we  assume  that  the  abundance  for  all  low-FIP  elements 
is  3  x  its  photospheric  value,  and  the  abundance  for  all  high-FIP 
elements  is  equal  to  its  photospheric  value  (Grevesse  et  al.  2007). 
For  the  LOS  depth,  we  adopt  the  value  of  0.05  RQ  at  1.1  RQ  so 
the  LOS  depth  at  a  given  height  r  is  0.05  x  r/ 1 . 1  Rq  under  spher¬ 
ical  expansion.  For  individual  CME  events,  this  LOS  depth  can 
be  estimated  from  the  length  and  inclination  of  the  neutral  line 
at  the  eruption  site  (Ciaravella  &  Raymond  2008;  Vrsnak  et  al. 
2009;  Lin  et  al.  2009).  In  any  case,  different  assumptions  for  the 


634 


KO  ET  AL. 


Vol.  722 


Height  (Rs) 


Height  (Rs) 


Height  (Rs) 


Height  (Rs) 


Figure  8.  Ionic  fraction  profiles  for  DR  height  at  1.2  R&  for  the  Streamline  model.  Plotted  are  for  three  coronal  models,  N1T1  (green  dash),  N1T2  (black  solid),  and 
N2T2  (red  dash  dot). 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 


Figure  9.  Electron  density  profiles  for  DR  height  at  1 . 1  (black)  and  1 .5  (red)  Rq 
for  the  Streamline  model.  Plotted  are  for  three  coronal  models,  N1T1  (dashed), 
N1T2  (solid),  and  N2T2  (dash-dotted).  Also  plotted  for  comparison  are  ne 
derived  from  UVCS  and  LASCO  data  (see  the  caption  of  Figure  4). 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 


FIP  bias  and  LOS  depth  can  easily  be  taken  into  account  since 
the  calculated  line  fluxes  are  linearly  proportional  to  these  two 
quantities. 

Figure  10  plots  the  predicted  fluxes  for  Sixn  A499,  [Caxiv] 
A943,  and  [Fexvm]  A.974  versus  heliocentric  heights  for  the 
Fully-mixed  model  with  three  coronal  models  and  DR  height 
at  1.2  Rq.  As  previously  mentioned,  these  lines  are  chosen 
because  they  were  commonly  observed  by  UVCS  in  post-CME 
CS  events.  UVCS  data  for  four  events  are  also  plotted  for 
comparison.  Figure  11  plots  the  same  but  for  the  Streamline 
model.  The  two  models  predict  similar  line  fluxes  at  low  heights 
(<  2.0  Rq).  The  difference  is  at  higher  heights  where  the  line 
fluxes  in  the  Fully-mixed  model  decrease  more  rapidly  with 
height.  At  5  RG,  these  lines  are  fainter  by  an  order  of  magnitude 
in  the  Fully-mixed  model  than  the  Streamline  model.  At  heights 
of  the  available  UVCS  observations  (^1.7  RQ),  the  model 
predictions  agree  with  the  observations  in  an  overall  sense.  The 
1998  March  23  event  seems  to  be  more  consistent  with  the 
N1T1/N1T2  models,  and  the  2002  January  8  event  seems  to 
be  more  consistent  with  the  N2T2  model.  This  would  indicate 
different  ambient  coronal  conditions  for  different  events.  For 
the  2003  November  4  event,  however,  Si  xn  and  Fe  xvm  lines 
are  consistent  with  different  coronal  models.  One  possible 
explanation  is  that  none  of  the  coronal  models  adopted  here 
represents  the  actual  coronal  condition  for  this  event.  A  higher 
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Figure  10.  For  the  Fully-mixed  model:  three  line  fluxes  vs.  height  for  three  coronal  models,  N1T1  (green  dash),  N1T2  (black  solid),  and  N2T2  (red  dash  dot)  for  DR 
height  at  1.2  Rq.  Also  plotted  are  the  data  from  UVCS  for  the  CS  events  on  1998  March  23  (orange)  at  1.50  Rq,  2002  January  8  (blue)  at  1.53  Rq,  2002  November 
26  (green;  shown  is  the  range  from  November  27,  00:26  UT  to  November  29,  00:20  UT)  at  1.61  Rq,  and  2003  November  4  (purple;  three  points  from  November  4, 
21:06  UT  to  November  5,  00:12  UT)  at  1.66  Rq.  The  coronal  abundance  is  assumed  to  be  3x  and  1  x  the  photospheric  value  for  low-FIP  and  high-FIP  elements, 
respectively,  and  the  LOS  depth  is  assumed  to  be  0.05  Rq  at  1.1  Rq.  Note  that  for  the  Sixn  flux,  that  of  the  2002  January  8  event  has  the  background  corona  taken 
out.  Others  should  be  regarded  as  upper  limit. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 
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Figure  11.  Same  as  Figure  10,  but  for  the  Streamline  model. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 


temperature  in  CS  (e.g.,  from  smaller  coronal  ,  Equation  (2)) 
would  produce  less  Sixn  emission  relative  to  the  Fexvm 
emission.  For  Ca  xiv,  all  models  underestimate  the  line  emission 
up  to  a  factor  of  5.  This  line  is  blended  with  [Si  vm]  A944  but 
the  absence  of  low-temperature  emission  along  the  slit  where 
the  Fe  xviii  emission  exists  rules  out  line  blending  as  a  possible 
cause.  The  discrepancy  may  lie  in  the  atomic  data  and  the 
assumed  coronal  abundance. 

Figure  12  shows  an  example  of  how  different  streamlines 
contribute  to  the  line  fluxes  observed  at  height  r  =  1 .7  Rq  for 
the  case  of  tdr  =  1.2  RQ.  All  three  coronal  models  are  plotted 
for  comparison.  Depending  on  the  individual  line,  the  dominant 
emission  comes  from  “streamlines”  of  a  limited  range  of  the 
SMS  crossings.  There  is  little  difference  among  the  three  coronal 
models.  This  is  in  contrast  with  the  ionic  fractions  (Figure  7), 
probably  because  of  a  combination  of  the  effects  of  the  electron 
density  and  temperature  on  the  line  emission. 

One  interesting  aspect  to  examine  is  the  change  of  line 
emission  with  time.  As  the  reconnection  continues,  the  location 
of  the  DR  will  move  up  toward  higher  height  with  speed 
approximately  equal  to  the  inflow  speed  of  the  reconnection 
(Lin  &  Forbes  2000).  Since  the  formation  of  the  SMS  is  much 
faster  (in  Alfven  speed)  than  the  movement  of  the  DR,  the  SMS 
developed  from  a  given  DR  height  can  be  regarded  as  a  steady- 
state  structure.  The  properties  within  the  CS  will  then  change 


with  the  movement  of  the  DR.  The  time  evolution  of  these 
properties  may  provide  information  about  the  DR  movement, 
therefore  the  inflow  speed/reconnection  rate.  Figures  13  and  14 
plot  the  three  line  fluxes  and  the  electron  density  versus  DR 
height  as  seen  at  three  locations  (1.5  RQ,  1.6  Rq,  1.7  Rq), 
along  with  the  UVCS  data,  for  the  Fully-mixed  model  and 
Streamline  model,  respectively.  The  time  sequence  of  line  fluxes 
and  electron  density  in  the  2003  November  4  event  (purple 
symbols;  Ciaravella  &  Raymond  2008;  observed  at  1.66  RQ) 
indicates  that  during  2.5  hr,  the  Fe  xviii  fluxes  decreased  by  a 
factor  of  0.63,  Sixn  fluxes  increased  by  a  factor  of  1.5,  and 
the  electron  density  decreased  by  a  factor  of  0.7.  If  we  use 
model  N1T2  as  a  proxy  (note  that  the  direction  of  change 
in  these  three  quantities  is  the  same  as  the  data,  even  though 
the  magnitude  does  not  all  fit),  this  roughly  corresponds  to  an 
upward  motion  of  around  0.15  RQ  which  in  turn  corresponds  to 
an  inflow  speed  of  ~10  km  s-1  and  MA  ~0.01  (based  on  the 
Alfven  speed  of  ~O(103)  km  s-1).  This  is  consistent  with  both 
theoretical  (although  on  the  low  side,  e.g.,  Lin  &  Forbes  2000) 
and  observational  (e.g.,  Yokoyama  et  al.  2002;  Ko  et  al.  2003; 
Lin  et  al.  2005)  findings. 

3.2.  Simulated  Broadband  X-ray  Emission 

Besides  particular  lines  of  interest  shown  above,  our  model  of 
the  electron  density,  temperature,  and  ionic  fractions  for  the  13 
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Figure  12.  Contribution  to  the  total  flux  observed  at  height  1 .7  Rq  from  streamlines  entering  the  SMS  at  heights  from  1 .2  Rq  (DR  height)  to  1.1  Rq.  Plotted  are  for 
three  lines  for  three  coronal  models,  N1T1  (green  dash),  N1T2  (black  solid),  and  N2T2  (red  dash  dot).  The  SMS  height  step  is  0.01  Rq. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 
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Figure  13.  For  the  Fully-mixed  model:  line  fluxes  and  electron  density  observed  at  1.5  Rq,  1.6  Rq,  1.7  Rq  vs.  the  height  of  the  diffusion  region  for  three  coronal 
models,  N1T1  (green  dash),  N1T2  (black  solid),  and  N2T2  (red  dash-dot).  Also  plotted  are  the  UVCS  data  (see  the  caption  of  Figure  10).  For  the  three  curves 
corresponding  to  the  three  observed  heights  of  each  model,  lower  height  has  higher  flux/density  in  general  (cf.  Figure  10). 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 


most  abundant  elements  enables  us  to  calculate  emission  spectra 
across  a  continuous  wavelength  range.  One  such  application 
is  to  use  the  emission  code  of  Raymond  &  Smith  (1977; 
updated  in  Cox  &  Raymond  1985)  which  calculates  the  line  and 
continuum  emission  spectrum.  This  predicted  spectrum  can  then 
be  compared  with  observations  from  broadband  instruments 
such  as  Hinode/X RT  and  SDO/AIA,  after  folding  with  the 


instrument’s  wavelength  response  function,  for  specific  CS 
events  these  instruments  might  observe. 

To  demonstrate  this,  we  use  the  non-equilibrium  ionization 
states  computed  with  both  models  along  with  the  elemental 
abundances  described  above  to  compute  the  emission  over  the 
range  2.5-4998  eV  (2.48-4960  A).  The  Raymond  &  Smith 
model  includes  fewer  emission  lines  than  the  APEC  code  (Smith 
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Figure  14.  Same  as  Figure  11,  but  for  the  Streamline  model. 
(A  color  version  of  this  figure  is  available  in  the  online  journal.) 


et  al.  2001),  but  for  low  spectral  resolution  data  the  results 
are  very  similar.  We  then  fold  the  emission  spectrum  with  the 
wavelength  response  function  of  the  XRT  instrument  on  Hinode 
to  obtain  the  predicted  count  rate  of  the  CS  if  it  were  to  be  seen 
by  XRT.  The  wavelength  response  function  is  calculated  using 
the  standard  SolarSoft  routines  with  contamination  thickness  for 
2008  April  9,  14  UT.  We  chose  this  particular  instrument  and 
date  because  a  very  interesting  CME/CS  event  was  observed  by 
SOHO  and  Hinode  on  this  day  that  has  been  analyzed  in  great 
detail  (Savage  et  al.  2010;  Landi  et  al.  2010). 

Figures  15  plots  the  predicted  count  rate  (DN  s-1  pixel-1) 
versus  height  for  the  Fully-mixed  model  as  seen  by  Hinode /XRT 
thin  Al-Poly  filter  for  the  three  coronal  models  and  three  DR 
heights.  Figure  16  plots  the  same  for  the  Streamline  model.  The 
LOS  depth  is  assumed  to  be  the  thickness  of  the  CS,  5  x  103  km, 
as  determined  by  Savage  et  al.  (2010)  for  the  2008  April  9  event. 
We  can  see  that  the  count  rate  can  differ  by  an  order  of  magnitude 
among  coronal  models  and  is  mainly  dictated  by  the  ambient 
coronal  density,  not  the  electron  temperature.  Figure  17  shows 
the  CS  observed  by  Hinode /X RT  on  2008  April  9.  Plotted  in 
the  upper  panels  are  the  measured  count  rates  along  the  solar  Im¬ 
position  at  four  solar  X-positions  that  cross  the  CS  (marked  in  the 
lower  panels).  The  location  of  the  DR  is  at  ~  1 .25  RQ  based  on 
the  downflows /upflows  measured  by  Savage  et  al.  (2010).  This 
is  around  the  second  outermost  X-position  marked  in  Figure  17 


(inner  solid  line).  The  “bumps”  between  Y pixel  of  100-300  are 
where  the  CS  is  located.  Plotted  in  red  on  the  two  lowest  curves 
(solid  lines)  are  the  Y  pixel  ranges  chosen  to  estimate  the  mean 
and  standard  deviation  of  the  count  rate  within  the  CS.  These 
numbers  are  shown  as  blue  data  points  in  Figures  15  and  16. 
The  heliocentric  heights  of  these  data  points  were  estimated 
based  on  that  the  flare/CME  source  was  ~  23°  behind  the  limb 
at  14  UT  on  2008  April  9  (Savage  et  al.  2010).  Comparing  the 
predicted  and  observed  count  rates  at  these  X-positions  indicates 
that  our  model  calculations  are  in  the  agreeable  range  with  the 
data  for  the  N1T1  and  N1T2  models,  but  not  for  the  N2T2  model. 
Note  that  our  predicted  count  rates,  bracketing  the  data  between 
DR  height  of  1.2  RQ  and  1.3  RQ,  are  also  consistent  with  the 
position  of  the  X  point  (~  1.25  RQ)  as  found  out  by  Savage 
et  al.  (2010).  The  agreement  implies  that  the  ambient  coronal 
density  profile  is  close  to  the  “Nl”  profile,  at  least  at  around 
1.3  Rq.  It  should  be  possible  to  compare  the  “Nl”  profile  with 
the  electron  density  profile  derived  from  the  WL  polarization 
brightness  measurement  at  the  time  before  and  during  the  event 
and  see  if  they  agree,  but  we  do  not  make  such  attempt  in  this 
paper. 

3.3.  Solar  Wind  Charge  States  at  1  AU 

Many  of  the  magnetic  clouds  (MCs)  seen  as  interplanetary 
coronal  mass  ejections  (ICMEs)  bear  high  ion  charge  states 
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Figure  15.  Predicted  count  rate  vs.  height  for  Hinode/XKT  thin  Al-Poly  filter 
for  the  Fully-mixed  model.  Plotted  are  for  three  coronal  models,  N1T1  (green 
dash),  N1T2  (black  solid),  and  N2T2  (red  dash  dot),  and  for  three  DR  heights 
at  1.1  Rq,  1.2  Rq,  1.3  Rq.  The  LOS  depth  is  assumed  to  be  5  x  103  km.  Also 
plotted  (blue  data  points)  are  four  measurements  of  Hinode/X RT  on  the  CS 
event  on  2008  April  9  (see  Figure  17). 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 

usually  represented  in  high  0+7/0+6  ratios  (>  0.2)  and  high- 
ionization  charge  states  of  Fe,  e.g.,  Fe+16  (Fepri  &  Zurbuchen 
2004).  Since  the  ion  charge  states  usually  freeze-in  near  the  Sun, 
this  implies  electron  heating  in  the  CME  material  in  the  early 
stage  of  the  eruption  process  (e.g.,  Akmal  et  al.  2001;  Rakowski 
et  al.  2007;  Fee  et  al.  2009).  Our  model  of  the  post-CME  CS  also 
predicts  higher  charge  states  due  to  higher  electron  temperature 
than  the  “normal”  corona  from  the  electron  heating  crossing  the 
SMS,  as  shown  in  our  results.  Plasma  ejected  upward  in  the 
CS  enters  the  flux  rope  formed  by  reconnection  below  the  CS 
core.  If  there  is  a  pre-existing  flux  rope,  reconnection  forms  a 
flux  rope  around  it  having  a  comparable  amount  of  magnetic 
flux,  and  therefore  volume  (Fin  et  al.  2004;  Mostl  et  al.  2008), 
so  that  much  of  the  plasma  in  a  MC  has  passed  through  the 
CS.  However,  there  is  no  report  so  far  of  definite  signatures  of 
post-CME  CS  measured  in  situ  at  1  AU.  Note  that  this  “CS” 
is  different  from  the  “in  situ  reconnection  exhaust”  reported  by 
Gosling  et  al.  (2005),  but  the  model  work  presented  here  can 
in  principle  be  applied  to  this  case  with  very  different  in  situ 
ambient  conditions  outside  the  SMS. 

Figures  18  and  19  show  the  predicted  frozen-in  0+7/0+6 
ratios  and  the  average  Fe  charge  at  1  AU  for  the  Fully-mixed 
model  and  Streamline  model,  respectively.  We  can  see  that  there 
are  notable  differences  between  the  two  model  approaches,  as 
well  as  the  three  coronal  models.  In  general,  the  Streamline 
model  predicts  higher  charge  states  at  1  AU  (or  >20  RQ 
according  to  Figures  6  and  8)  than  the  Fully-mixed  model. 
The  ionization  charge  state  is  higher  for  model  N2T2  than 
the  other  two  coronal  models  due  to  higher  CS  temperature 
(e.g.,  Figure  4),  as  opposed  to  that  the  emission  is  lower  for 
model  N2T2  due  to  lower  density.  We  note  that  the  range  of 
these  predicted  quantities  of  0+7/0+6  and  average  Fe  charge  is 
common  in  the  ICME  solar  wind  data  (0+7/0+6  >  0.2,  average 
Fe  charge  >12;  see,  e.g.,  Fepri  &  Zurbuchen  2004,  Wimmer- 
Schweingruber  et  al.  2006,  and  Zurbuchen  &  Richardson  2006). 


Figure  16.  Predicted  count  rate  vs.  height  for  Hinode/XKT  thin  Al-Poly  filter 
for  the  Streamline  model.  See  the  caption  of  Figure  15. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 

If  some  of  the  ICMEs  with  high  0+7/0+6  and  Fe  charge  states 
are  associated  with  such  post-CME  CS  structure,  it  would  imply 
that  the  DR  remains  low  in  the  corona.  It  remains  to  be  seen 
if  the  models  presented  here  can  be  tested  by  observations  of 
charge  states  within  the  to-be-identified  post-CME  CS  structure 
in  the  solar  wind. 

4.  DISCUSSION 

In  order  to  understand  post-CME  WF  rays  and  CSs,  it  is 
important  to  have  quantitative  comparisons  between  theory 
and  observations.  The  models  presented  here  are  a  first  step 
in  that  direction.  Several  crucial  observations  exist,  including 
the  Fe  xxiv  spines  seen  above  some  post-flare  loops  in  TRACE 
images  (Innes  et  al.  2003a),  the  high  velocity  Fe  xix  and  Fe  xxi 
lines  seen  by  SUMER  in  several  events  (Innes  et  al.  2003b; 
Wang  et  al.  2007),  the  UVCS  observations  of  highly  ionized 
gas  (Ciaravella  et  al.  2002;  Ko  et  al.  2003;  Bemporad  et  al. 
2006;  Ciaravella  &  Raymond  2008),  evidence  of  reconnection 
downflows/upflows  observed  by  Hinode/XKT  (Savage  et  al. 
2010),  and  in  situ  measurements  of  high  ionization  states  in 
ICMEs  (e.g.,  Fepri  &  Zurbuchen  2004). 

We  have  taken  a  simple  model  of  the  exhaust  region  of 
Petschek  reconnection  and  computed  the  time-dependent  ion¬ 
ization  state  in  order  to  predict  UV  emission  line  intensities 
and  the  count  rates  in  X-ray  bands.  This  model  is  consistent 
with  observations,  thus  it  is  a  viable  explanation  for  the  ob¬ 
served  “post-CME  current  sheet”.  This  does  not  exclude  other 
possible  models  but  any  model  should  be  able  to  predict  prop¬ 
erties  that  can  produce  observables  that  are  consistent  with  the 
observations.  The  Vrsnak  et  al.  (2009)  model,  in  our  opinion, 
is  the  most  viable  model  so  far  that  adequately  treats  plasma 
properties  of  the  CS  above  the  DR  for  which  we  can  calculate 
physical  parameters  and  compare  with  the  observations.  Given 
the  model  assumptions,  the  good  agreement  with  the  observa¬ 
tions  in  both  density  and  UV  line/X-ray  emission  is  more  than 
just  fortuitous.  In  the  context  of  this  model,  the  width  of  the  WF 
ray  increases  with  height  which  is  also  consistent  with  the  data 
(see  Vrsnak  et  al.  2009).  This  indicates  that  this  model  can  be 
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Figure  17.  Lower  panels:  Hinode/XRT  Al-Poly  image  on  2008  April  9  at  13:46  UT  (lower  left)  and  17:06  UT  (lower  right)  showing  the  CS  extended  above  the 
post-CME  loops/cusp.  The  vertical  lines  are  the  solar  X-positions  chosen  for  the  plots  in  the  upper  panels.  The  spatial  binning  is  1  pixel.  Upper  panels:  observed  count 
rate  in  DN  s-1  along  the  F-pixel  number  at  the  four  solar  X-positions  marked  by  the  vertical  lines  in  the  XRT  images.  The  “bumps”  between  Y  pixel  of  100-300  are 
where  the  CS  is  located.  Note  that,  according  to  Savage  et  al.  (2010),  the  DR  is  located  in  between  the  dashed  and  solid  lines.  Thus  our  model  results,  which  only 
apply  to  the  CS/SMS  above  the  DR,  should  be  compared  only  with  the  data  at  the  locations  marked  by  the  two  solid  lines  (lowest  two  curves  in  the  upper  panels). 
Plotted  in  red  on  the  two  lowest  curves  are  the  T-pixel  ranges  used  to  estimate  the  mean  and  standard  deviation  of  the  count  rate  within  the  CS  (shown  as  blue  data 
points  in  Figures  15  and  16). 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 


the  right  interpretation  of  the  observed  WL  ray/CS.  In  principle, 
our  calculations  can  be  applied  to  any  post-CME  CS  models  that 
provide  information  about  the  electron  density /temperature  and 
ion  flow  speed,  such  as  in  the  MHD  models  (e.g.,  Riley  et  al. 
2008  by  using  proton  parameters).  Note  that  the  main  purpose 
of  this  paper  is  to  show  how  this  particular  post-CME  CS  model 
concept  agrees  with  the  observations  but  we  do  not  attempt  to 
make  fine  adjustments  of  the  model  parameters  to  fit  data  of 
any  one  particular  CS  event.  It  is  also  possible  to  apply  such 
calculations  to  model  the  downflow  region  below  the  DR  (e.g., 
Yokoyama  &  Shibata  1998,  2001)  but  we  do  not  make  such 
attempt  in  this  paper. 

Both  CS  models  (“Fully-mixed”  and  “Streamline”)  are  found 
to  give  similar  UV  and  X-ray  emission  at  heights  below  ~2  RQ 
but  it  drops  more  quickly  with  height  in  the  Fully-mixed 
model.  Thus,  it  seems  that  CS  observations  higher  than  2  RQ 


would  be  able  to  differentiate  the  two  models.  However,  current 
instruments  do  not  have  the  sensitivity  yet  to  detect  such  high- 
temperature  emission  much  beyond  what  is  presented  here. 
Figures  4  and  9  indicate  that  electron  densities  derived  from  WL 
observations  beyond  4  RQ  can  be  used  to  differentiate  the  two 
models.  Such  data  are  available  although  we  do  not  attempt  to 
acquire  in  this  work.  The  ionic  charge  state  (Figures  1 8  and  19)  is 
another  useful  parameter  to  differentiate  these  models.  However, 
an  unambiguous,  positive  identification  of  such  post-CME  CS 
in  the  solar  wind  is  necessary  before  such  comparison  can  be 
made. 

The  model  results  depend  strongly  on  the  model  of  the 
ambient  corona.  The  differences  among  the  three  coronal 
models  show  that,  within  the  context  of  the  models  here,  the 
coronal  electron  density  profile  is  the  most  important  factor 
in  determining  the  emission  properties  inside  the  CS.  In  this 
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Figure  18.  0+7/0+6  ratio  and  average  charge  of  Fe  observed  at  1  AU  for  the 
Fully-mixed  model  vs.  height  of  the  DR  with  the  three  coronal  models  compared. 
(A  color  version  of  this  figure  is  available  in  the  online  journal.) 


Height  of  Diffusion  Region  (Rs) 


Height  of  Diffusion  Region  (Rs) 
Figure  19.  Same  as  Figure  18,  but  for  the  Streamline  model. 
(A  color  version  of  this  figure  is  available  in  the  online  journal.) 


work,  we  did  not  investigate  the  effect  of  the  ambient  coronal 
magnetic  field.  But  we  expect  it  will  also  play  a  crucial  role  (as 
dictating  the  plasma  /3)  in  determining  the  CS  properties.  The 
more  accurately  the  ambient  conditions  at  the  time  of  the  CS 
event  are  known,  the  better  such  observations  can  pinpoint  the 
reconnection  processes  and  properties  in  the  post-CME  CS.  And 
comparison  of  such  models  presented  here  with  observations  can 
offer  a  possible  diagnostic  tool  for  inferring  the  external  coronal 
conditions  of  the  post-CME  CS. 

As  a  final  remark,  the  calculation  of  plasma  emission  and 
time-dependent  ionization  relies  on  atomic  rates  for  collisional 
excitation,  ionization,  and  recombination.  Aside  from  the  ac¬ 
curacy  issue  of  these  atomic  rates,  the  electron  velocity  dis¬ 
tribution  function  (VDF)  is  an  important  factor  in  these  rate 
calculations.  In  this  work,  we  assume  that  the  electron  VDF  is 
a  Maxwellian.  The  deviation  of  the  VDF  from  a  Maxwellian 
can  affect  these  rates  in  different  degrees  for  different  ions  (e.g., 
Ko  et  al.  1996;  Dzifcakova  &  Mason  2008).  One  future  effort 
of  our  calculations  will  be  to  investigate  possible  effects  of  a 
non-Mawellian  electron  VDF.  A  kinetic  approach  in  magnetic 
reconnection/CS  models  that  can  provide  information  about  the 
electron  VDF  would  be  helpful  in  nailing  down  our  knowledge 
of  the  post-CME  CS  with  observations. 
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